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Abstract 

We have measured the absolute branching ratio of the tt'^tt^ (7) decay, using 

~20 milhon tagged mesons cohected with the KLOE detector at DA<I>NE, the 
Frascati (/>-factory. Signal counts are obtained from the fit of the distribution of the 
momentum of the charged decay particle in the kaon rest frame. The result, inclusive 
of final-state radiation, is BR(K+ 7r+7r° (7))=0.2065 ± O.OOOSstat ± O.OOOSsyst- 

Key words: e+e- Experiments, Kaon decays 
PACS: 13.25.Es 



1 Introduction 

The branching ratio of the — > tt+tt" (7) decay (-^'7^2) is pai't of the KLOE 
program of precise and fully inclusive kaon branching ratios (BRs) measure- 
ment. We have already measured the main Kl pll2] and Ks [3P] branching 
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ratios. We report here our measurement of the BIi{K~^ tt+tt'' (7)) which 
together with BR{K^ [5], BR{K^ ttH^u) [6j and BR(K± ^ 

Tr^Tr^vr") [7] covers 95% of all charged kaon decays. The importance of this mea- 
surement is twofold: i) the most recent measurement based on 16,000 events 
from a sample of ~10^ kaon decays, BR{K+ 7r+7r° (7))=0.2118 ± 0.0028 [8] , 
dates back to more than 30 years ago and gives no information on the radi- 
ation cut-off and ii) this BR is necessary to obtain BR(i^'i3) from measure- 
ments normalized to BR(i^'^2) PfTO] - The — »• 7r"'"7r° (7) branching ratio 
can be used, together with the Ks — *■ tttt branching ratios [3], to determine 
the relative phase 6q — 62 of the 1=0 and 1=2 s-wave Trvr-scattering ampli- 
tudes [H] . In the following we report our measurement of the absolute branch- 
ing ratio BR(i^"'" — ^ tt+tt" (7)) performed with the KLOE detector using an 
integrated luminosity /Ldt~250 pb~^ collected at DA$NE, the Frascati 0- 
factory. DA$NE is an e~^e~ collider operated at the energy of 1020 MeV, the 
mass of the meson. Equal energy positron and electron beams collide at 
an angle of (tt— 0.025) radians producing 0-mesons with a transverse momen- 
tum of ~13 MeV. In its rest frame, the 0-meson decays into anti-collinear 
K^K~ pairs of ~127 MeV momentum and this remains approximately true 
in the laboratory. Detection of a K"^ (the tagging kaon) therefore signals the 
presence of a (the tagged kaon) of given momentum and direction. This 
procedure, called tagging, allows measurements of absolute BRs. 



2 The KLOE detector 

The KLOE detector consists of a large volume drift chamber surrounded by 
an electromagnetic sampling calorimeter. The entire detector is immersed in 
an axial magnetic field B = 0.52 T. The drift chamber (DC) [12], 3.3 m long 
and 4 m in diameter, has a stereo geometry with 12,582 drift cells arranged 
in 58 layers and operates with a 90% helium-10% isobutane gas mixture. 
Tracking in the DC provides measurements of the momentum of charged par- 
ticles with cr(p_L)/p_L < 0.4% for polar angles larger than 45°. The spatial 
resolution is ~150 /^m in the bending plane, ~2 mm on the z coordinate 
and ~3 mm on decay vertices. The electromagnetic calorimeter (EMC) ^3] 
consists of a cylindrical barrel and two endcaps, covering a solid angle of 
98% of 47r. Particles crossing the lead-scintillator-fiber structure of the EMC, 
segmented into five planes in depth, are detected as local energy deposits. 
Deposits close in time and space are grouped into clusters. The energy and 
time resolution for electromagnetic showers are oeIE = ^.7Vo/ \J E{GeV) and 

at = 57 ^ E{GeV) © 100 ps, respectively. The trigger [14j requires two 
isolated energy deposits in the EMC with: E > 50 MeV in the barrel and 
E > 150 MeV in the endcaps. Cosmic-ray muons are identified as events with 
two energy deposits with E > 30 MeV in the outermost EMC planes and ve- 
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toed at the trigger level (CRV). A software filter (SF), based on the topology 
and multiplicity of EMC clusters and DC hits, is applied to reject machine 
background. The effect of both CRV and SF on the BR measurement must be 
determined. In the following the coordinate system is defined with the 2;-axis 
along the bisector of the e^e~ beams, the y-axis vertical and the x-axis toward 
the center of the collider rings and origin at the collision point. 



3 The measurement 



Tagging with K" ji'v {K~2) and K~ 7r"7r° {K~2) decays provides 
two samples of pure for signal search. These two-body decays are easily 
identified as peaks in the distribution of the p* variable, the momentum of the 
charged decay particle in the kaon rest frame evaluated using the pion mass, as 
described in |T5j. The tagging kaon is required to satisfy the trigger request by 
itself, minimizing the dependence of the trigger efficiency on the decay mode 
of the tagged kaon. The residual dependency, which we refer to as the tag 
bias in the following, must be determined for the BR evaluation. We choose to 
measure BIi{K.,^2) using mesons because for them the nuclear interaction 
correction is negligible, since the probability of interaction is ~10~^ for 
and -3.4% for R-. 

The branching ratio is determined as: 



BR (A- ^ W) = X \ , (1) 

where Nk+^-k+ttO {-y) is the signal count, A^Tag the number of tagged events 
and e is the overall efficiency, including the detector acceptance edet and the 
reconstruction efficiency erec- The detector acceptance (edet~59%), entering in 
the final efficiency evaluation, is taken from MC and its value is related to the 
charged kaon lifetime r. Consequently the BR depends on r as: 

BR(") /BR(°) = 1 - 0.0395 ns-^(r - r(°)) (2) 



with T^^^ = 12.385 ns, the current world average value [IS]. A variation of the 
lifetime of 0.1% changes the BR of 0.05% of its value. The corrections Cqry, 
CsF and Ctb account for the cosmic-ray muons veto, the software filter and 
tag bias effects, respectively. 

The sample used for this measurement has been processed and filtered with 
the KLOE standard reconstruction software and event classification proce- 
dure [IS]. The KLOE Monte Carlo (MC) simulation package, GEANFI, has 
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Fig. 1. Fit of the distribution. Left: black dots are data and grey histogram is 
the fit output. Right: the three contributions used to fit the data are shown: i^^2j 
^7^2 and three-body decays. 

been used to produce an event sample equivalent to the data. The different 
operating conditions of DA$NE during data taking, machine parameters and 
background, are included in the MC on a run-by-run basis. The simulation also 
includes final-state radiation [T^ guaranteeing correct measurement of fully 
inclusive BRs. The result of the simulation has been compared with theoretical 
predictions and experimental results for several kaon decay channels. 

3.1 K^2't(^99^d sample 

The number of tagged by K~2 decays, the /r~2-tagged sample, is A^xag = 
12,113,686. The signal selection uses DC information only. The K'^ track 
is identified as a positive track with point of closest approach (PCA) to the 

interaction point (IP) satisfying Y^^|cA"+~yfcA < 10 cm and |-2pca| < 20 cm, 
and momentum 70 < < 130 MeV. The PCA is evaluated extrapolating the 
track backwards to the IP taking into account energy losses. Decay vertices 
(V) are accepted in the fiducial volume 40 < ^Jx'^ + y'^ < 150 cm, |zv| <150 
cm. Loose cuts on p* and on the difference between the momenta of the kaon 
and the charged secondary track, 50 < p* < 370 MeV and —320 < Ap < —50 
MeV, reject K 3tt decays. 

The K^2 signal count is extracted from the fit of the p* distribution (Fig. [1]). 
The p* spectrum has two peaks: the first at ~236 MeV due to muons from 
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K^2 decays, and the second at ~205 MeV due to pions from Kt,2 decays. The 
contribution from three-body decays shows at lower p* values. Having used 
the pion mass for the evaluation, the K^2 peak is distorted. We fit the p* 
distribution between 180 and 350 MeV using three contributions: -^7r2 and 
three-body decays. The shapes of the K^2 and K.„2 peaks are obtained from 
data control samples, selected using EMC information only. The i^^2 spectrum 
is obtained from the i^'^2-control-sample used for the efficiency evaluation and 
described later. The Kf^2 spectrum is obtained from the control sample selected 
for the BR(i^'"'" fJ'~^^) measurement |5J. Once a tagging K~2 decay has been 
identified, we ask for only one EMC cluster with energy -Eciu > 80 MeV and 
no clusters with energy between 20 and 80 MeV. There are no requirements 
on EMC clusters with energy below 20 MeV, in order to retain — /i"""// 7 
decays and ii^v decays with machine background clusters in the EMC. 

This high-purity sample (~99%) is called ii'^2-control-sample. Bin by bin MC 
corrections account for small distortions induced in the p* spectra by the 
control sample selections. The three-body component is obtained from the 
MC simulation, which has been tuned with the data K^2- and i('7r2-control- 
sample. Fig.[T]left shows the result of the fit of the p* distribution compared to 
the data, while the three different contributions are visible on the right. The 
fit, with xVndf = 154.5/156 (P(x^) =0.52), gives iVx+^^+7rO (^) = 818, 347 ± 
1,912, the error accounting for the statistics (not only data). 

The reconstruction efficiency tree has been evaluated with data. Since the K^2 
events are identified from DC information, the data control sample is selected 
using EMC information. Once a tagging K~2 decay has been identified, we 
construct by kinematics the track from the K~ track. We then search for 
two photons in the EMC and, using their time and energy information, we 
determine the decay point, the di-photon mass and momentum. The best 
accuracy is obtained minimizing the sum of the square of the differences be- 
tween the decay time from photons and path and between the di-photon 
mass and the vr" mass. Having determined the track and decay point of the 
and the vr*^ direction, we determine the expected vr"*" track using the two- 
body decay hypothesis. The kinematics of this hypothesis is then verified by 
requiring the presence of a cluster in the EMC, with a distance from the pion 
track (iciu < 30 cm. These events define the i^7r2-control-sample. The con- 
tamination from decays without a 7r° in the final state is ~ 0.1%. About 
5% contamination from Ki^^ decays is present and becomes about 3% after 
signal selection. Corrections accounting for small distortions due to the selec- 
tion of the data control sample have been evaluated using MC. Defining etrue 
the true efficiency to reconstruct signal decays in the DC volume and the 
reconstruction efficiency obtained using the i^',r2-control-sample, the average 
correction to be applied to the efficiency is etmc/ecs ~ 0.99. The efficiency to 
be used in eq. HJis e = 0.3176 ± 0.0005. 

The corrections Ccmv = 1.0005 ±0.0003 and C^p = 1.0183 ± 0.0003 have been 
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htatisticai errors 


— " J. 

using K^^ tag 


Statistical errors 


using tag 


Source 


Value (%) 


Source 


Value (%) 


fit signal count 


0.23 


fit signal count 


0.27 


efficiency 


0.12 


efficiency 


0.13 


SF and CRV 


0.04 


SF and CRV 


0.03 


e correction 


0.12 


e correction 


0.13 


TB 


0.05 


TB 


0.05 


Total 


0.30 


Total 


0.33 



Table 1 

Summary of fractional statistical uncertainties on BR(ir^2) measured using K^2' 
and i^^2"tagged samples. Left K~2 tag, right K~2 tag. 

measured with data taken without the cosmic-ray muons veto and the software 
filter, respectively. The correction for the tag bias, Ctb = 1.0106 ± O.OOOSstat, 
has been evaluated using MC. The distributions of variables used for the 
selection of the tagging decay have been checked with data. Table [1] left lists 
the statistical fractional uncertainties on the branching ratio measurement and 
the total value is 0.3%. 

3.2 K~2-t(igged sample 

The number of tagged by K~2 decays, the i^'~2"tagged sample, is iVxag = 
9, 352, 915. Table [2] compares the values of the Cqrw, Csf and Ctb corrections 
obtained for K~2- and i^'~2"tagged events. The two tags have very different cor- 





Kf,2 tag 


K^2 tag 


CcRV 


1.0005±0.0003 


1.0007±0.0003 


Csf 


1.0183±0.0003 


1.00093±0.00006 


Ctb 


1.0106±0.0005 


1.009±0.0006 



Table 2 

Corrections to BR(^C^2) measured using K~2- and -R'~2"tagged samples. 

rections for the effect of the software filter (SF). The Csf correction measured 
using the K~2 tag is ~1.8% while using the K~2 tag is ~0.1%. The same signal 
selection as before is applied to the sample tagged by K^2 decays and the fit 
of the distribution determines the signal count. The spectra used for the fit 
have been obtained as described in the previous section, once a tagging K~2 
decay has been identified. The signal count is iV^+^^+^o = 621, 612 ±1, 678. 
For the efficiency evaluation we have used the K^2-control-sample tagged by 
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O 4- J- ■ 

Systematic errors 


using tag 


C J- 4-' 

bystematic errors 


usmg tag 


hource 


Value (%) 


hource 


Value (%) 


p* fit range 


0.06 


fit range 


0.07 


shape 


0.12 


Kfj_2 shape 


0.14 


^7^2 shape 


0.16 


Kt^2 shape 


0.17 


efficiency 


0.30 


efficiency 


0.30 


rv 


0.17 


rv 


0.17 


lifetime r 


0.12 


lifetime r 


0.12 


TB 


0.01 


TB 


0.01 


Nucl. int. 


< 0.02 


Nucl. int. 


< 0.02 


Total 


0.42 


Total 


0.43 



Table 3 

Summary of fractional systematic uncertainties on BR(i^^2) measured using K^2' 
and i^^2"tagged samples. Left K~2 tag, right tag. 

K~2 decays. The efficiency is e = 0.3182±0.0005, corrected for the control sam- 
ple selection and the detector acceptance taken from MC. The total statistical 
fractional uncertainty on BR(i^'~2) measured using the i^'~2-tagged sample is 
0.33%. Table [T] right summarizes the fractional statistical uncertainties. 



4 Systematic uncertainties 

The systematic uncertainties on BR{K^2) from K~2 and /r~2"tagged samples 
are listed in tables [3] left and right, respectively. The stability of both BR 
measurements with respect to different data taking periods and conditions has 
been checked. A detailed discussion of the systematic studies follows. These 
studies have been done varying the selection cuts in wide intervals and checking 
the stability of the BR. 

The lower bound of the p* fit range, 180 MeV, has been moved from 165 to 
190 MeV, changing by almost a factor of two the contribution from three- 
body decays. For each value of the lower bound, we have performed the fit of 
the distribution and evaluated the overall efficiency. We observe a minimal 
change of the BR value in the above range. The maximum variation of the 
BR is taken as systematic uncertainty. The contributions to the fractional 
systematic uncertainty on the BR are 0.06% {K~2 tag) and 0.07% {K~2 tag). 

The spectrum of the K^2 component for the fit of the p* distribution is ob- 



8 



tained from the i^^^2-control-sample. The Kf^2 spectrum is most affected by 
the cut at 20 MeV on the cluster energy Eqiu [IS], connected to the accep- 
tance of a photon from — > /i+t/ 7 decays or from machine background 
events. The stabihty of the BR measurement has been checked by changing 
the -Eciu cut from 10 to 30 MeV, corresponding to a change in the purity of 
the i^'^2-control-sample from ~99.3% to ~97%. Neghgible effects are observed 
with -Eciu values larger than 30 MeV. The maximum variation of the BR has 
been taken as systematic uncertainty. The fractional systematic uncertainties 
are 0.12% {K'^ tag) and 0.14% {K'^ tag). 

The spectrum of the component for the fit of the distribution is ob- 
tained from the i^7r2-control-sample. The systematic effect has been estimated 
performing the fit with the p* spectrum obtained from a different control 
sample. We select decays in the DC, using the signal selection of sec. 13. Ij 
and require the identification of a 7r°, looking for two photons in the EMC 
fulfilling the following requests. The two photons have to be on-time: the dif- 
ference between the kaon decay times, evaluated using the cluster time and 
the distance between the decay vertex and the cluster position, has to be 
within 3(Tt (see sec. The kaon decay time from the kaon path and from the 
photons have to be compatible within resolutions. The difference between the 
di-photon mass and the 7r° mass has to be within 3cr, with (T~18 MeV. The 
Ki3 contamination of this sample is 20%, larger than the 3% contamination of 
the i^7r2-control-sample. Thus the spectrum of the Kt,2 component obtained 
from this sample needs larger MC bin by bin corrections (as large as 60%) 
compared to the default used (20% at maximum and for low values of p*). 
Using this spectrum we have been performed the fit of the p* distribution, 
also varying the fit range. The BR results are in agreement, within errors, 
with the values obtained using the spectrum from the i^'7r2-control-sample. 
The maximum difference between the BRs obtained with the two spectra has 
been taken as systematic uncertainty. The fractional contribution is 0.16% 
(ir;2 tag) and 0.17% (R-^ tag). 

The reconstruction efficiency has been evaluated with data using the 
control-sample. The systematic uncertainty has been estimated using a differ- 
ent control sample, with larger Ki^ contamination (~11% compared to ~3%) 
and MC correction to be applied to the efficiency (~12% compared to ~1%). 

decays with a vr*^ in the final state are selected, as done for the ii',r2-control- 
sample but without the dc\u cut. We determine the p* of the charged secondary 
track, using the two-body hypothesis and the and vr*^ momenta. Two-body 
decays are then selected applying the asymmetric cut 0.5a < p* — 205 < a, 
with cr~18 MeV, around the peak at 205 MeV of the p* distribution. The 
BRs measured using the efficiencies obtained from the above sample and the 
i^7r2-control-sample agree within errors. Conservatively the difference between 
these two BRs is taken as the systematic uncertainty. The contribution to the 
fractional systematic uncertainty is 0.3%. 
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The decay vertex has to satisfy the requirement 40 < pv = y a^v + 2/v < 150 
cm. The lower bound of the Pv range, p™™ = 40 cm, has been moved from 38 
to 42 cm with the detector acceptance changing of ~6% of its value. For each 
p™" value, we have performed the fit, evaluated the efficiency and measured 
the BR. The efficiency has been evaluated with the ii'^2-control-sample. The 
resolution on the decay point, using only the time information in the EMC, 
is cr~1.5 cm. Thus the above interval corresponds to a change of more than 2cr. 
The BR results are in agreement within the statistical error and their rms is 
taken as systematic uncertainty. The contribution to the fractional systematic 
uncertainty is 0.17%. 

The BR depends on the charged kaon lifetime r through the detector accep- 
tance. The systematic effect has been obtained using eq. [2] and the 0.24% 
fractional accuracy of the KLOE measurement r=12.347±0.030 ns [20j. The 
contribution to the fractional systematic uncertainty is 0.12%. 

The fractional systematic uncertainty from the tag definition is 0.01%, as 
obtained changing separately the requirements to identify the tagging decay. 

The fraction of undergoing nuclear interaction has been evaluated using 
the MC simulation and considered as upper bound value of the systematic un- 
certainty. The contribution to the fractional systematic uncertainty is <0.02%. 

The emission of radiation in the 7r+7r°7 decay is dominated by the 

Inner Bremsstrahlung (IB) contribution and the Direct Emission term can be 
neglected [2T]. The predicted value of the IB branching ratio is 2.61x10"^ 
and our simulation reproduces this value within few 10~^. Thus, the sys- 
tematic uncertainty associated to the modelling of final-state radiation is 
negligible. The average of the experimental results is BR(i^"'" — >■ 7r"'"7r°7) = 
(2.75±0.15) X 10"'^ [TS] and its error gives a negligible contribution to the un- 
certainty on BR(_ft'^2)- 

The fit of the p* distribution providing the count for K^2 decays, gives the 
number of — > p+i/ (7) decays as well. The reliability of the fit procedure is 
confirmed by comparing BR(ir^2) finding agreement with our published 
result [5]. The criteria for signal selection and efficiency evaluation from this 
reference have been followed. There is therefore a correlation p{K^2, ^^^2) of 
—3.4% between our BR(i^'^2) BR(ii'^2) measurements using the signal 
count extracted from the fit procedure. For our published BR{K^2) result [5] 
we did not use the fit of the p* distribution to extract the signal count. The 
number of ^ fi^u (7) decays was obtained by counting the number of 
events with p* > 225 MeV, after background subtraction. Therefore there is 
no correlation between the published BR(i^^2) value and our BR(_ft'^2) value. 

When averaging the BR(ii'^2) values obtained from the K~2- and i^'^2"t^gg6d 
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samples we have to account for correlations. The same data control sample 
for efficiency evaluation has been used for both measurements, thus giving 
a correlation in the statistical as well as in the systematic contribution to 
the BR uncertainty. The contribution to the systematic uncertainty from the 
charged kaon lifetime value r is common to both measurements as well as the 
contribution from the value. The correlation between the two BK{K^2) 
measurement is 56%. 



5 Conclusions 



We have measured the branching ratio of the —>■ tt+tt^ (7) decay, fully 
inclusive of final-state radiation, using samples tagged by K~2 and K~2 
decays. From 12,113,686 ii'~2"tagged events, we find N^+^^+^o = 818, 347 ± 
1, 912 signal counts. Using eq. [T] we obtain the branching ratio: 

BR 7r+7r°(7))|^ = 0.20638 ± 0.00062stat ± 0.00087syst- (3) 



From 9,352,915 i^^2-tagged events we have Nk+- 
signal counts corresponding to: 



621,612 ± 1,678 



BR {K^ 



TT 



-7r°(7)) 



K^2-tag 



0.20668 ± 0.00068,t«t ± 0.00089, 



^syst • 



(4) 



The above BRs are evaluated using the current average value for the 
lifetime r(°) = 12.385 ns (see eqj2]). Averaging these two results, accounting 
for correlations, we obtain: 

BR(ir+ ^ 7r+7r°(7)) = 0.2065 ± 0.0005stat ± 0.0008,ysf (5) 



This absolute branching ratio measurement is fully inclusive of final-state ra- 
diation and has a 0.46% accuracy. Our result is 1.3% (~2o") lower than the 
PDG fit [in]. The global fit to all available charged kaon measurements [TS] 
gives BR(i^'"'" — > 7r"'"7r''(7)) = 0.2064 ± 0.0008, in agreement with our result. 

We fit the six largest BRs and the lifetime r using our measurements of 
r [20], BR(ir+2) (eq-ED, BR(ir+2) 13, ^^{K^) [6J and BR{K^ n^n^^) [7], 
with their dependence on r, together with BR(i^'^ — > tt^tt+tt") from the 
PDG04 averagj^ [22J , with the sum of the BRs constrained to unity. The fit 
results, with x^/ndf = 0.59/1 (CL=44%), are shown in tablelHand confirm the 
validity of our measurement (eq. E]), assuming the correctness of BR(i^'^ 
TT^'^Tr+Tr^). 

^ PDG '06 gives the result of their constrained fit but not the average of the data 
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Parameter Value 



Correlation coefficients 



0.6376(12) 
0.2071(9) 



+0.48 



BR(7r±7r+7r-) 0.0553(9) 
BR{Ke3) 0.0498(5) 
BR{K^3) 0.0324(4) 



+0.34 



+0.37 



-0.48 



+0.21 



-0.12 



-0.13 



+0.16 



+0.15 +0.58 



BR(7r=^7r°7r°) 0.01765(25) -0.11 
r (ns) 12.344(29) -0.15 



+0.05 



-0.21 



-0.07 -0.06 -0.05 -0.015 



0.05 +0.04 +0.04 



Table 4 

Results of the fit to BRs. 

We can also evaluate BR^K"^ 7r=''7r+7r^) by using our measurements of the 
above listed BRs and imposing the constraint J2BR{K^ ^ /) = 1. With 
BR(i^'+2) = 0.63660 ± 0.00175, BR{K^ tt^tt^tt^) = 0.01763 ± 0.00025 and 
BR{K^2) = 0.20681 ± 0.00094, BRiK^^) = 0.04972 ± 0.00053 and BR{K%) = 
0.03237 ± 0.00039, evaluated at r equal to our measured value 12.347±0.030 
ns, we get BR(i^'^ tt^tt^tt^) = 0.0568 ± 0.0022. This result is in agreement 
with the PDG04 average BR{K^ ^ vr^vr+Tr-) = 0.0550 ± 0.0010 [52]. 

The K —>■ data can be used to extract the s-wave vrvr scattering phase 
shift difference 6o — 62 a-t s = m|^. This has been done by the FlaviaNet Kaon 
Working Group [23] , using as experimental inputs the FlaviaNet world average 
value of the BR for tt+tt^ (7) decay [IB] and the KLOE measurement 

of the BRs for Ks tt^tt^ and Ks n~^7r~{'y) decays [^. Their result is 
So — ^2 = (57.5 ± 3.4)°, accounting for strong and electromagnetic isospin 
breaking effects. The same calculation has been done using our value of the 
BR for K+ 7r+7r° (7) decay from the fit to K"^ BRs (table S]). The result 
is Sq — 62 = (57.6 ± 3.4)° and confirms the ~2(T discrepancy with the results 
from phenomenological analysis of tttc scattering amplitudes [23J. 
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